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Abstract 

We study an extension of one dimensional Calogero model involving strongly coupled and elec- 
trically charged particles. Besides Calogero term ^f^-, there is an extra factor described by a Yukawa 
like coupling modeling short distance interactions. Mimicking Calogero analysis and using develop- 
ments in formal series of the wave function ^ (x) factorised as (x) with e{e — 1) — g, we develop 
a technique to approach the spectrum of the generalized system and show that information on full 
spectrum is captured by $ (x) and $" {x) at the singular point a; = of the potential. Convergence 
of f dx\'^ {^)\^ requires e > — ^ and is shown to be sensitive to the zero mode of $ {x) at a; = 0. 
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1 Introduction 

The study of integrability of hamiltonian systems has been one of the basic tasks during decades 
as it gives precise informations on classical and quantum dynamics as well as spectrum of the hamil- 
tonians This analysis has lead to the discovery of several features regarding quantum spectrum 
of integrable hamiltonian systems The familiar examples are given by the remarkable class of 

integrable low dimensional field models based on both finite dimensional and affine Lie algebras. The 
study of integrable hamiltonian systems has also shown that there is only a few number of real physical 
systems that are completely solvable. The usual harmonic oscillator and one dimensional Calogero model 
as well as extensions certainly some of these examples. Despite the variety of methods used 

to study integrability, one may, roughly speaking, say that the main difficulty in getting exact solutions 
comes essentially from the fact that the underlying differential equations are non linear coupled relations 
difficult to solve. Lessons learnt from the study of hamiltonian systems and infinite dimensional sym- 
metries of integrable models let understand that hamiltonian systems may in general be classified into 
three basic sets: (i) integrable hamiltonian systems where the spectrum of the hamiltonian is exactly 
determined as in case of Calogero model IH], (ii) quasi-integrable systems where, though not completely 
determined, the hamiltonian spectrum, obtained by perturbation of exact models, is under control and 
(iii) remaining others with unknown spectrum. 

Motivated by remarkable properties of integrable one dimensional Calogero type interaction , 
we study in this paper the case where Calogero particles are moreover electrically charged and strongly 
coupled. Coulomb interaction |_^/';^ | and strong coupling between pairs of particles {xi,Xj} are imple- 
mented by adjunction to the usual Calogero hamiltonian an extra short distance interaction described 
by a Yukawa like coupling exp (—a \xi — Xj\) / \xi — Xj \ where a is the inverse of Debye length A. Among 
our results, we show that for the case of two particles, the discrete spectrum {En, \E'n (x) = a;'^$„ (x)} of 
this generalized Calogero system with e > — | is completely determined by the knowledge of the wave 
factor <i>„ (x) and $^ {x) at the singular point a; = of the potential. The range e > — i is required by 
the convergence of / dx \'^n and is shown to be sensitive to the zero mode of $„ (x) at x = 0. 

The organisation of this paper is as follows: In section 2, we review briefly quantum Calogero 
model by solving directly the Calogero wave equation for the case of two particles; the case of N particles 
is straightforward and follows the same lines. Comments regarding this solution are also given; they are 
used in the interpretation of our generalization. In section 3, we study the extension of Calogero model 
involving charged and strongly correlated particles described by a Yukawa like interaction. Motivated 
by similarities with Calogero analysis and inspired from moment technique, we develop a method to 
approach the discrete spectrum of the generalized Calogero system. In section 4, we give our conclusion 
and make a discussion regarding some specific limits. 

2 ID Calogero model: brief review 

One dimensional Calogero model describing the quantum dynamics of two interacting particles, 
of local coordinates Xi and X2 and relative position x — {x2 — Xi) /v^, is governed by the following 
eigenvalue wave equation H'i' {x) — E'^ (x), 



where to fix the idea we take a; > 0; i.e X2 > xi. In this relation, uj is the usual harmonic oscillator 
frequency and g is the Calogero coupling parameter. To get the discrete spectrum {En, (x)} of this 




(1) 



differential equation, it is interesting to use the following direct method which is also valid for the case 
of a system of N particles. This method relies on solving directly eqQ) and involves three steps each one 
associated with the factorisation of the field variable \l/ (x) as a;^<i> (x) = x^P {x) exp (— ^a;^) for some 
modulus e. First, set the change (x) = x^<J> (cc), which for e > 0, reflects just the fact that two classical 
particles {xi,X2} cannot live together at xi — X2] the wave function should then vanish; i.e ^ (0) = 0; 
but $ (0) ^ 0. Note that the choice e > is not a necessary condition, since quantum mechanically we 
should rather have, 

/•OO 

dx\^ {x)\^ <oo. (2) 

As far convergence of this integral is concerned, one should not that there are two sources of possi- 
ble divergences, one for a; ^ and the other for a; ^ oo. Decomposing eq® as dx\^ {x)f' + 
dx \ ^ {x)'f' + dx {x)\^ for some small positive number y and a large positive z] then using the 
fact that /°° dx \^ {x)\^ is finite^, we will prove later that in the region < a; < y, convergence requires 
that 2e+i 1*^ C*^)! should be finite. By substitution of the change 5* (x) = a;'^(f> {x) in eqlU), we get the 
following constraint relation on $, 

$" (x) + — $' {x) + {2E - u'^x') $ {x) = 0, (3) 

X 

which, for a matter of discussion of the solution, we prefer to rewrite it as follows, 

{x) + —^' {x) + 2{E-U)<^{x) = Q, (4) 

X 

with U {x) = Lo^x'^/2. Note that the factorisation 5" (x) ~ x^$(a;) has been also introduced to fix 
the power e in term of Calogero coupling constant g 'as g — e {e — 1). This relation follows from the 
requirement that the change 5* = a;^$ has to kill the problematic term which is one of the two 
sources of difficulty (non linearity) in the search for exact solution. The other difficulty is treated in the 
second step which consist to make the local scaling change, 

$(a;) -P(a;)exp(-|x2) , (5) 

which for matter of generality, we rewrite it also as <i> (x) = P (x) exp (— / (x)) with / (x) — ^x^ = ^. 
Using this change, we can bring ea H3l4|l to the form, 

P" (x) + - 2/'^ P' (x) + (2E -2U- 2^/' - /" + P (x) = (6) 

By substituting 2U = /'^ = w^x^, we discover the following familiar differential equation, 

P" (x) + ( — - 2wx ) P' (x) + [2E~u; {2e + 1)] P (x) = 0, (7) 



where now the coefficient [2E — uj {2e + 1)] in front of P (x) is constant. This is a typical equation 
|15j whose solution is given by a special hypergeometric function which we describe in the next step. 
Setting u = tux'^ and P (x) = L (u) and putting back into above equation, we get the following familiar 
differential relation, 

(u) + - ^) in) + (I; - ^) L (u) = 0. (8) 

This differential equation has a known exact solution given by Laguerre polynomials L„ ^ (m) provided 
that the factor - 2£±i) is quantized, i.e - 2£±i) ^neN, 

[."+"e-1, (9) 



^For large x, say z < x < oo, Yukawa like potential may be neglected with respect to Calogero interaction. So 
/^°° dx I* (x)P ^ dx I {x)\^ which is known to be finite. 



where the upper index a is a convention notation; it should not be confused with a power. Therefore 
the discrete spectrum {En, (x)} of quantum Calogero two particle system reads as follows: 

*„(.t) = x'expi^-'^x^^Lir^^Lux^). (10) 

The ground state is given by ^Pq (x) — exp (— f 2:^) and i?™' = (e + 5) w; for e = (g = 0); one 
recovers the usual spectrum of the ground state for harmonic oscillator. 

For later use, it is also interesting to note the three following: (a) First remark that the above 
solution may be also rewritten in a positive modes series as follows, 

00 

®(-')(x;6,..)=^C^2t(e,c.) (11) 

fc=0 

where the expansion coefficients cl^2k ^) depending on the coupling moduli, 

^n,2k {^^^) - [jk^])M2^TuT7T^) ' ^^^^ 

with = n\/j\ [n — j)\ and where for later use note that C^^q''^ 7^ 0. Though forbidden by the reality 
condition of total wave ^ {x), the manifest discrete symmetry x <-> (—a;) of solution (|ll|l is not hazardous; 
it is a symmetry of the Calogero hamiltonian operator Under the exchange of the position xi ^ X2, 
the complex wave function maps as, 

^{xi,X2)=e'^'m{x2,xi), (13) 

showing that, for e different from even integer, the two particles obey a non trivial statistics, (b) The 
second thing we want to note is that the coefficients C^™'^ and C^™'^of the series ((TT|) are related as 
follows, 

^">2 -(27TT)^"'" ■ ^ ^ 

From ea Hll(l . one sees that above relation may be also rewritten as Ecai = — (2e + 1) ^^"^l^j . (c) The 
third thing we want to comment is that quantity (/'^ — 2^/' — /" — 2J7) appearing in eq® is constant, 
say equal to K. Its solution is obviously given by C/ = ^a;^, / = ^x^ and if = — (2e + 1). Note also if 
we take / = (e — 1) In a;, one gets precisely U = 1^ which is nothing but the first factorisation we have 
used to absorb Calogero term in the hamiltonian. 



3 Calogero model with a Yukawa like potential 

In this section, we consider the generalization of previous analysis to the case where the Calogero 
particles are electrically charged and are moreover strongly correlated. The interaction describing this 
particular behavior is given by the following Yukawa like potential |16| . 

Vrufe (a;) = ^ exp (^-^^ , a:V2 = |a:i - a;2| , (15) 

where a is a coupling constant and A is the Debye length. Note that for the range a; < A, the two leading 
terms of Vyuk {x) are respectively given by the usual Coulombian term ^ and a constant — The 
presence of the constant shows that energy spectrum E of this generalized system should read as 



E = Ecai ~ ^+ other contributions coming from Coulomb term and remaining polynomial fluctuations 
(exp (— f ) ^ 1 ~ f + ^ (f )^)- -'^'^^ a; > A, the spectrum of the quantum system is mainly given by the 
usual Calogero one exp (— f ) ^ 0). The Shrodinger equation describing the quantum dynamics of 
these interacting particles reads, in local coordinates xi and X2, as, 

^-^^^ + ^(a:i,X2)^ ^{xi,X2)=E^>{xi,X2), (16) 

where ^' is the total wave function, V is the total translation invariant potential given by, 

V(.) = -.^ + ^ + -^exp(^-_j, (17) 

and E is the total energy eigenvalue depending on the potential moduli uj, g, a and A. Note that, like 
for Calogero system, transformation generated by space scale dilatation x (^x, maps the hamiltonian 
Ti. as follows Ti. — > C, '^'H provided that couplings are mapped as w — > C^^^i 9 ^ 9i ^ C^^o: and 
A CA. Note moreover that the limit a small is also equivalent to large values of the relative distance 
X with respect to the Debye length A (x >> A). In this region, the additional Yukawa like interaction 
exp (^—^x^ may be neglected. 

To get the discrete spectrum {En, ^'n (x)} of the wave equation eas (|16ll7|l . we shall follow the 
Calogero method outlined in previous section and the moment technique to be introduced later on. 
Using classical properties of the wave function and translation invariance (xi — > + est), we can thus 
decompose ^ {x) in three factors as before 5* (x) = x^ exp {—^x'^) P {x) where a priori e is positive; but 
quantum effects requires that it should be as e > — ^ as we show later, see also eq(j21l and J^l- In the first 
step, the factorisation of 5* reads as follows ^' [x) — x'^®{x). Like in Calogero analysis, this factorisation 
is also used to kill the non linear term ^ by making an appropriate choice of e parameter {e — e (.9)). 
Finding 4" (x) is then equivalent to determining the e parameter in terms of potential moduli and the 
unknown function ®{x). Substituting 5-" = x"®" + 2ea;'^"^®' + e(e - l)^^^^® in eqlO, Schrodinger 
equation becomes. 



2e , 

— e + 

X 



e(£-l)-g 2 2 2V2a 
Lo X exp 



(^) 


+ 2E 







= 0. (18) 



This is a second order differential equation with non constant coefficients which, as we know, its solution 
isn't a simple matter. However, ea (jl8|l may be simplified a little bit by working, like in Calogero model, 
in a special region of the coupling constant moduli space; 

e(e — 1) — g = 0, a, A free moduli. (19) 

Note that g — g (e) is quadratic function in e and so it has two zeros namely e_ {g = 0) = and 
£+((7 = 0) = 1, see also figure, As e is the basic parameter in this analysis, it is then interesting to 

express it like e = e (g); but this leads however to two possible solutions e± (g) given by. 



1- yiq^ 1 + 

e-(ff) = 2 ' ^+(3) = 2 ■ ^ ' 

Generally speaking, one should distinguish between these two cases; however this not necessary since 
these values are related as e+ (g) + e_ {g) = 1. One derives the quantum spectrum for one of the cases, 
say e = e+ (5), then deduces the other one by using the identity e_ = 1 — e+. The second step is to 




Figure 1: This figure describes the variation of Calogero coupling constant g in term of the statistical parameter 
e. As one sees there is two regions for Calogero coupling g, attractive and repulsive. In the interface g — , we 
have two solutions (£-,<?) = (0,0) and (£+,(?) = (1,0), which, by ea ll'^ . shows that they correspond to Bose 
and Fermi statistics respectively. For generic values of g, we have the two branches e± (g) = iiV^+S^ 



factor $ as $ (x) = P {x) exp {—^x'^) in order to absorb the harmonic term cj^x^. By substitution, the 
previous equation becomes 

P" +(-- 2xu;) P'+(2E~i2e+l)cj~ ^^e"^ ) P = 0, (21) 



where obviously e stands for e± (g) — i^^^^Ss, Strictly speaking we have two wave equations to solve 
and then two discrete spectrum to get; we shall refer to them as {4*^, E^}. Note moreover that 
the two above factorisation of the field variable 5* {x) implies, amongst others, the two following: (i) 
The first change of variables fixes the parameter e as e± = whose reality property requires 

that g > — J showing that there are various phases, an attractive one associated to the moduli space 
— i < ,9 < and which may be interpreted as a confining phase, two repulsive ones given by g > 
with — i < e_ < and e+ > 1 and finally two interfaces associated with the points (e-,5) = (0,0) 
and (e+,g) — (1,0). Note that quantum mechanically, the classical argument 5* (a;) =0 for a; = 
leading to e > is not a quantum principle. The right condition is to require rather that 5* (x) is a 
Hilbert space wave function with a convergent integral J dx j^* (x)|^ < 00. Following |lfi| . wave functions 
such as 5* (x) — x^<i> (x) , with e > are also Hilbert space solutions despite the apparent divergence 
at X = 0. This feature is a typical quantum mechanical property which may be derived in several 
manners; one of which is the elaborated one developed in using canonical operator formalism; a 
second one is given by the direct method to be developed below and leading to ea lt^ . To have a 
quick idea on the rule behind this e shift to left, we present here below a naive way to see it. The 
argument relies on considering ea (|21|l and requires vanishing of fundamental energy for arbitrary values 
of relative position x. This is a typical feature dictated by renormalisability of quantum systems and 
allow to redefine energy origin by absorbing the constant term. Recall that in quantum mechanics, the 
fundamental energy constant term captures quantum effects and is one of the sources of divergences. 
In canonical quantization for instance, this term has some thing to do with operators normal ordering 
and is absorbed in the renormalisation process shifting energy as i? > 0. For our concern, it is enough 
to look at boundaries of the system (|x| — > 00) where Yukawa like interaction may be ignored. The 
corresponding energy positivity condition _B > is then the same as in Calogero case and reads as 
^ > ^£+1 equivalently 2e + 1 > 0. In general, one can prove that for the case of N particles, this 
condition extends as Ne + 1 > showing that we should have e > ^ . In what follow we shall give an 



other way to deal with this quantum shift, (ii) The second factorisation of the wave function tells us 
that as far as moduli space of coupling constants is concerned, one may also consider the moduli space 
region where the term exp (^^x^) ^^"^ treated as a perturbation around Calogero potential 
(j << 1). There, the spectrum {E^ (e, w, a. A) ,Vn (2;; e, w, a. A)} we are looking for should reduce to 
Calogero one, plus fluctuations as shown below; 



V, 



(x;e,c.,a,A) - pi--') ^a:;e,Lu) + O (^a,j 

i?„(e,a,A) = Ei''^''>{e,u;) + o(^a,jy (22) 

where Vn°'^\x]e^uj) and En"'^^ are the exact Calogero spectrum ea (|10f) : see also conclusion section. 
The normalisation condition of the total wave function ^P, which, by using the factorisation 5* (x) = 
x^P(x) exp(— ^x^), reads as, 

1 



Jo 



'f{x)\Ux^ I x^'P^ {x)e-'^'' dx, e>-2' ^^^^ 



requires convergence of the integral. For x = 00, such condition is ensured provided that x^^P^ {x) < 
e~"^ and for x = 0, the quantity ^2e+i (0) should be finite. Both these conditions are fulfilled for 
Laguerre polynomials associated with Calogero model. In the case where Yukawa like coupling are 
implemented, the function P {x) is no longer polynomial but has expansion type ^g^^PsX'^ involving 
positive modes only; that is with a good behavior near x = 0. Under this assumption we can prove 
that ea (|23|l converges as well for the range ^ > —\ exactly as in Calogero system. Let us give some 
details on the proof. First note there are two sources for possible divergences of ea (|23|) . one at x = +00 
and the other at x = 0. At infinity, Yukawa like potential ^e~"^ is negligible and the system is mainly 
the Calogero one where, as noted before, solutions are known for the range e > — i. In this case P (x) 
is a polynomial eq(ll) and its convergence is ensured by the gaussian factor whatever the value of e 
is. At X = however, the situation is some how subtle; we have e~'^^ ^ 1 — ujx^ + O (x'') and so a 
possible divergence may come from low modes of the expansion of P (x) which may be approximated 
in this region as po + Pix + O (x^) independently on whether P {x) is a polynom or not; the only 
requirement is that P (x) should have positive modes only. Therefore, we have x^'^P^ (x) e~'^^^ ~ 
PqX^' + 2poPix'^'^^^ + O (x^^+^), which diverges for e < ^ since, 

dx C ^2_^2.+l ^Pom 2.+2 PI-PO^ 2.+3 ^ q / (24) 

' ^ ^' 2e + l 2e + 2 2e + 3 ^ ' ^ ' 

At 2e + 1 = 0, one has a logarithmic divergence; it corresponds to a critical point describing collapsing 
particles as interpreted in For 2e + 2 = 0, one has both linear and logarithmic divergences. However 
in the range e > -j- the situation is healthy like in Calogero model. Note that convergence of/|«'(x)rdx 
is also sensitive to the values of the modes Ps- For po — for instance, even the logarithmic divergence 
at e = ^ disappear and, as one sees on ea lf^ . is shifted to the moduli space point e = 

Theorem 1 ; 

(a) The discrete energy eigenvalues En of generalized Calogero system with electrically charged particle 
and strongly correlated described, in addition to Calogero term, by a Yukawa like interaction eq^^ are 
given by 

£;„(e,a,A) = -(26+l)|i|| + /(a,A) (25) 

where $„ (0) and <I>^ (0) are the values of the wave function $ (x) and its second derivative at x = and 
where f (a. A) is equal to H£(^£:i£l^ 



(b) The total wave functions (a;) — x'^^n (x) of this generalized Calogero system is completely deter- 
mined up to the knowledge of the field (zero mode) values {x — 0) and $" (x = 0). 

Having in mind the method of solving standard Calogero wave equation presented in section 
2 and convergence of H23|) for P {x) with an expansion type X]s>o^'s^^ involving positive modes, one 
may proceed ahead to establish this theorem by mimicking the same method. The idea is to start from 
ea H18|l which we rewrite as, 



2e 



2En - uj'^x'^ 



2V20 



■ exp 



'X^/2 



= 0, 



P" {x) +{^- '2^^^ P' {x) + {2E - (2e + 1) w - C/) P (x) 



(26) 



where U — e ^ and try to put it in a form type (|7I8|I . 

To prove the theorem, we shall develop a method inspired from the solution and an expansion 
of the wave function as a formal series. The key idea may be summarized as follows: first use the 
feature da; |\I>„ (a;)|^ < 00 ea (l23|l and the factorisation (a;) = a;^P„ exp (— ^a;^) to build the 
wave functions ^'„. Second benefit from known properties of Calogero wave function, in particular its 
development as a series (|11I12II with positive modes which we suppose also valid in presence of Yukawa 
like interaction, 

*»(a;) = ^C„,.(e,w,a,A) a;^ e > 0, (27) 

r>0 

where the Cn,k coefficients, depending on the coupling moduli, are the new unknown quantities which 
have to be determined. This means that, like vE'^^"'-' (a;) = x*^''^"'^ (x), the wave (a;) — x'$„ (x) may 
be expanded in a formal series as. 



(a;) = ^ Cn^r (e, Lo, a. A) 



(28) 



r>0 



To that purpose, note that in a distribution sense, the integral /j^^ :^^^x"^dx — J^-^ d (a;™+^) is usually 
equal to zero except for the case m = — 1 where the integral is different from zero. Using this property, 
which may be thought of as a real version of Cauchy theorem of complex analysis, we can inverse ea (|28f) 
as follows, 



Cn,i (e,w, a. A) 



<f„ (a;) = (0) , ;>0. 



(29) 



To get the explicit expression of the new variables Cn^i (e, w, a, A), we start from the differential equation, 

-xV^^ 



x^®';, + 2e {x®'J + 



2E„x'^ - up-x'^ 



aa;exp 



A 



0. 



(30) 



Putting the expansion ®„ (x) = X^fcLo ^n,k ^'^ back into ea H3U|) and using the development exp ^ ^x^ ^ i 



namely X]m=o ^ \"^m\ after term rearranging the following algebraic equation. 







(31) 



/£>0 



where the first four An.k coefficients are as Anfi = [e(e — 1) — 5] Cn.o, An,i = (2e + e(e — 1) — (/) Cn,i 
2\/2aCn.o together with 



Ar^ 



Ar. 



(2 (2e + 1) + e(e -\)-g) C„,2 - 2 
(3 (2e + 2) + e(e - 1) - g) C„,3 - 2 



- \ En 



2a 
T 

2a 
T 



n,0 



(32) 



^/2c 



and the remaining others (k > 4) read as — [k (2e + fc — 1) + e(e — 1) — g] Cn,k — Dn,k, with Dn,kS 
given by, 

-Dn.fe = 2a > Cn.fc-m-i + LU Cn.k-4 — '^EnCn.k-2 (33) 

m=0 

The solution of ea H31|) for arbitrary x is obtained by requiring the vanishing of the A„ j. coefficients 
which give in turn the following recurrent constraint eqs on the C„ ^ functions, 

= [e(e-l)-5]C„,o, 
_ 2aV2 
" 2e + e(e-l)-.9^"'°' 

2[V2aC„,i-(i;„ + ^)a,o] 
' 2{2e+l) + eie-l)-g ' (^'^ 



n,3 



2\/2aC„,2 - 2 + ^) C„,i + 2^C„,o 
3(e + l) + e(e-l)-5 



and for fc > 4, 



0„,fc = ■— — ; — — — [6b) 



1+1 

^ ' " ^ '2EnCn,k-2 

k{2e + k-l) + e{e-l)- g 
From these recurrent relations, one sees that in the case e 7^ the elements of the infinite set {Cn,k} are 
all of them function of Cn,o, the zero mode of the expansion H28|l . and the energy En- Notice that in the 
case e = 0, one is left with Yukawa like model since classically there is no Calogero interaction. In this 
situation and using the relations An.o — [e(e ~ 1) ^ ff] C'n,0: ^n,i = (2e + e(e — 1) — g) C„.i — 2\/2aC„^o 
together with An,2 and An,3 given by eqslj23), have either Cn,o ^ 0; but a = or Cn,o ~ and 
a =^ 0. In the second case, one should have Cn,i ^ otherwise the wave function vanishes. Non trivial 
solution^ requires Cn,o 7^ 0; the first equation of (|34|l fixes then the Calogero g coupling constant as in 
the case of absence of Yukawa like coupling, 

e{e - 1) = g. (36) 
Putting this relation back into above eas H34l35|) . we get the following reduced constraint eqs, 

C„,i 



aV2 
e 






2e + 1 V 


e A 




2a' 


3(e + l) 


e(2e + 



C„,2 = :^[—-=^-Ej, (37) 

_ (3e+l)(i;„ + ^) 1 

- , , e(2e + l) + 

together with the following recurrent relations, 

- /fc(2e + fc-l) ■ ^^^> 

From eas (|37|l . one sees that for e 7^ one can express the energy En in terms of the coefficients Cn,o 
and Cn,2 as shown below, 

f2a^ 2a\ , , C„ 2 

The knowledge of C„_o and C„.2 fixes then completely the spectrum {^/^ (x) , En}; but we have no exact 
way to determine C„_o and Cn^, except the use of approximation methods. This is the case where the 

■^The moduli space poles at 2t = (1 — fc), fc = 2, 3, should be compared with those of eQ l24i . 



correlation between particles are strong enough {x << A) so that the Yukawa like coupling may be 
interpreted as a perturbation with respect to Calogero interaction. In this limit, the ratio may be 

approximated by ea (fT^ and consequently the approximated energy reads as £'„ — + E^"'' 

with e ^ 0. As one sees Ai?„ = AE"™' ^ AiJ°'"^ whatever positive integer n is; it would be interesting 
to check this behaviour by other methods. Progress in this matter will be exposed in a future occasion. 
Note finally that the case e = with non trivial wave function requires Cn,o = 0, but C'n,i 7^ and the 
previous energy relation get replaced by 2a'^ ~ ^ ~ 

4 Conclusion and discussions 

In this paper, we have studied an extension of one dimensional Calogero model involving electrically 
charged particles with strong correlations. These couplings are described, in addition to Calogero 
potential by a Yukawa like interaction Vyuk (x) = ^ exp (— f ) , = \xi — X2\. Though this 
system is not completely integrable; we have shown that its discrete spectrum {'i'n (x) ~ a;^$„ (x) , En}, 
e > may be totally determined up to fixing the value of the wave factor $„ (x) and $^ (x) at the 
singularity a; = of the potential; that is the zero mode Cn,o and the second one C„.2 of the wave 
function expansion ea H28|l . The quantum condition e > ^ is shown to be intimately related with the 
non zero value $„ (0); see ea (l24ll . In the case of Calogero model, the modes C^^""}^, C™2 a-nd energy i?^"' 
are all of them fixed by a remarkable property of confluent hypergeometric functions which are solution 
of eq® and which, for discrete spectrum, reduces to Laguerre polynomials. In this study, we do not 
have a similar property due to the presence of highly non linear potential interactions (exp (— f ))■ We 
have focused on the case of two identical particles {xi, 2:2}; but this is just a matter of simplicity. Our 
analysis extends also to the case of a generic number N of identical particles {xi, xn}- There, the wave 
function ^E* depends on the ^^-^^~^'> relative positions = Xi—xj (5* (x) = {HijXijY $ {xi2, x^^i^n)) 
and the expansion of the factor $ {xi2, x^^i^pf) follows the same lines as for the case of two particles 
considered here. 

In the end, we would like to add two comments. First note that for g = e(e — 1), one should 
distinguish the moduli space regions: (i) — i < e < 0, (ii) < e < 1 and (iii) e > 1. Moreover 
solving g — e{e — l) as e± (g) eas (l^ . one sees that the corresponding spectrum {\E'^,i?^} may be 
split into {^2fe+i' ^2fc+i} ^'^^ {^2*:' ^2*:} with k positive integer. This splitting has a particularly 
interesting interpretation in the moduli space region g and a go to zero. Physically speaking, this 
corresponds to large x where the two potential components ^ and j^e a are negligible with respect 
to harmonic interaction ^x'^. In this limit, the discrete spectrum {'i'^,E^} should then tend to the 
usual spectrum of one dimensional harmonic oscillator. Setting g = a — Oin ea HlU|) . one 

gets En{e = 0,a = 0) = -E^"' (e = 0) = (2n + i) cj, n e N; which at first sight seems not reproducing 
the full energy spectrum of the harmonic oscillator that we recall here En'^'^^ = (n + i) w, n G N. The 
naive relation between E™' (e — 0) and En'^'^^ reads as, 

E^^'ie^O) = Eir\ (40) 

and shows that apparently there are missing terms associated with the odd integer modes E2°^^{- In 
fact there is no contradiction here since we have not implemented the full picture for the limit g = 0. In 

this case, one should distinguish two sectors €+ = 1,5 = and e_ = 0, g = whose respective energies 
(.9 = 0) and E,i"'"°'^ {g = 0) are given by hm^^o [(2n + 1) + ix/TT4^] oj = [(2n +l) + ^]uj 
and linig^o [(2n + 1) — ^-\/l + Ag] lu = [2n + -i] uj. They read in terms of En^'^'' as follows, 

(g = 0) = e'-I Eir-^^^^ (5 = 0)= EiT\ (41) 



As one sees Calogero interaction splits harmonic oscillator energies into two blocks £'2^+1 ^i'n'^^ 
with a global separation -E'2n+i ~ ~ ^ whatever the integer n is. This is a remarkable feature 

which should deserve more investigation. The second comment concerns the low mode structure of wave 
functions obtained from eas (|28l I!MI35|I by implementation of eas (|2()|l . In the limit e+i = \, g — Q 
and the case C^q ^ 0, we have for leading terms C,|^ = a^/2C^Q, C^2 = —j^{2a'^ ~ X ~ )' 
C+3 = ^(^(2a2 _ ^ - E+) - V2a{E+ + ^) + ^) together with the following modes = 
I(ITT) (2«E!rio ^'Zl^ Cl,_^_,+u^Cl,_,^2E+Cl,_,) for fc > 4. In the case = 0, 5 - 0; but 
= and + 0, we have = a^/2C-i, = ^(2a2 - f - ^^^,4 = %i(4^(2a2 - 
^ - E-) V2a{E- + ^) + ^) and for k > 4, C", = -^(Sa^^^i^ + 
'^'^^n fc-4 ^ "^^n^n k-2)- Notc that the modcs in the two cases are related as C^^(e+ = l)C^j(e_ = 
0) = C^j,_(^j(e_ — 0)C^o(^+ = last note that near x — 0, the wave function (x) = x'^^n (x) 

behaves as x'^ (^Cn,o + Cn,ix + O (x^)) . In the case where the leading modes are as Cn,o — and 
Cn,i 0, all happens as if e is shifted as ? = e + 1 and Cn,Q = Cn,i- 
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